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Abstract. We consider the differential equations y" = Xo(x)y' + so(x)y, where 
Xo(x),sq(x) are C°°— functions. We prove (i) if the differential equation, has a 
polynomial solution of degree n > 0, then 5 n = A rl s„_i — A rl _is„ = 0, where 

A„ = A^_ : + s n _i + A A„_i and s„ = s' n _ 1 + s X k -i, n = 1,2, Conversely 

(ii) if A„A n _i and 6 n = 0, then the differential equation has a polynomial solution 
of degree at most n. We show that the classical differential equations of Laguerre, 
Hermite, Legendre, Jacobi, Chebyshev (first and second kind), Gegenbauer, and the 
Hypergeometric type, etc, obey this criterion. Further, we find the polynomial solutions 
for the generalized Hermite, Laguerre, Legendre and Chebyshev differential equations. 
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1. Introduction 

The question as to whether a second order linear homogeneous differential equation 
has a polynomial solution has attracted much interest since the early classification of 
Bochner of othogonal polynomials In 1929, Bochner posed a problem of determining 
all families of orthogonal polynomials that are solutions of the differential equation 

(ax 2 + bx + c)y"(x) + (alx + e)y' n (x) - [i n y n [x) = 0. (1) 

Bochner found that, up to a linear change of variable, only the classical polynomials 
of Jacobi, Laguerre and Hermite and the Bessel polynomials satisfied a second-order 
differential equation [H] of the form (0). In general, the question as to which second 
order linear homogeneous differential equation has polynomial solutions (not necessary 
a sequence of orthogonal polynomials) is not easily answered, since it would involve 
studying a wide variety of equations, including those with regular and irregular singular 
points. In this article, we present a simple criterion for the existence of polynomial 
solutions of a differential equation of the form 

y" = X y' + s y (2) 

where A ,s are C°°— functions. A key feature of the present work is to note the 
invariant structure of the right-hand side of (J2J) under further differentiation. Indeed, if 
we differentiate (j2J) with respect to x, we find that 

y»> = \ iy > + Sl y ( 3 ) 

where Ai = Aq + sq + Aq and si = s' + soAo- If we find the second derivative of equation 
(JHJ), we obtain 

yW=\ 2 y f + s 2 y (4) 

where A 2 = X[ + s% + A Ai and s 2 = s[ + s Ai. Thus, for (n + l) th and (n + 2) th 
derivative, n = 1,2,..., we have 

y^ = K-iy' + sn-iy (5) 

and 

y(n + 2) = Xny > + Sny (6) 

respectively, where 

A n = A^_ x + + A A n _i and s n = s' n _ x + s Q A n _i. (7) 
From (jSJ) and © we have 

Ky {n+l) - K-iy {n+2) = S n y where 5 n = A„s n _i - A n _is„. (8) 

In an earlier paper ^T] we proved the principal theorem of the Asymptotic Iteration 
Method (AIM), namely 
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Theorem 1: Given X and Sq in C°°(a,b), the differential equation (HJ) has the general 
solution 

y(x) = exp I- J a{t)dt\ C 2 + C x J exp I j (A (r) + 2a(r))rfrj dt (0) 



a. (10) 



i//or some n > 

^ri S n —\ 

A n A n _i 

The present article is not about a classification of orthogonal polynomials which is 
well-established problem in the literature P0-[E]- Rather, the principal goal of the present 
paper is to characterize when Eq.(j2J) has a polynomial solution. In the next section we 
shall show that the differential equation (j2J) has a polynomial solution of degree n, if 
for some n > 0, 5 n = 0. In Section 3, we show through detailed analysis that the 
classical differential equations of Laguerre, Hermite, Legendre, Jacobi, Chebyshev (first 
and second kind), Gegenbauer, and the Hypergeometric type, etc., obey this criterion. 
In Section 4, we apply the criterion to obtain polynomial solutions to the generalized 
Hermite, Laguerre, Legendre and Chebyshev differential equations. As we shall show, 
the criterion presented here works whether or not the differential equation (j2J) has a set 
of orthogonal polynomial solutions, or a class of orthogonal polynomial solutions in the 
quasi-definite sense [TU| . 



2. A criterion for polynomial solutions 

The existence of polynomial solutions is characterized by the vanishing of 5 n . This is 
the principal theoretical result of this paper. We have: 

Theorem 2: (i) If the second-order differential equation ^) has a polynomial solution 
of degree n, then 

AnSn-l — A„_is„ = 8 n = 0. (11) 

Conversely (ii) if A n A n _x ^ 0, and 5 n = 0, then the differential equation (0) has a 
polynomial solution whose degree is at most n. 

Proof: (i) For the given differential equation (j2J), if y is a polynomial of degree at 
most n we have y( n+1 *> = y( n+2 ) = 0. consequently we conclude from (JHJ) that 5 n = 0. 
(ii) Conversely, if 5 n = and A„A n _i ^ 0, then we have s n _i/A n _i = s n /\ n = a, and, 

X 

from Theorem 1, we conclude that a solution is given by y = exp(— J a{t)di). Therefore, 
in particular, y' = —ay = —t^U- Consequently, from ?/ n+1 ) = A n _iw' + s n _iw, we infer 

A n _i 

that |/( Tl+1 ) = 0, or, equivalently, that y is a polynomial of degree at most n. □ 

Theorem 2 gives us the condition under which the given differential equation has 
a polynomial solution. Theorem (1), in particular ©, provides a tool for the explicit 
computation of these polynomials. In the next section, we apply these results to a 
variety of classes of differential equation: in each case we provide the explicit condition 
which yields polynomial solutions. 
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In this section, we apply Theorem 2 to the classical differential equations of 
mathematical physics. First, we give an alternative proof to Bochner's results (JTJ, 
using the criterion developed in Theorem 2. 

Theorem 3: The second-order differential equation (QJj has a polynomial solution of 
degree n if 

fi n = n(d + (n-l)a), n = 0,1,2,.... (12) 
The corresponding polynomial solutions are 

2/o=l 

iji = dx + e 

y 2 = (d + a)(d + 2a)x 2 + 2(6 + e)(d + a)x + e(b + e) + c(d + 2a) 
y 3 = (d + 2a) (d + 3a) (d + 4a)x 3 + 3(d + 2a) (d + 3a) (e + 2b)x 2 

+ 3(d + 2a)(6(3e + 26) + c(4a + d) + e 2 )x 

+ 4d6c + e 3 + 3dec + lOaec + 2e6 2 + 3e 2 6 



Proof: By means of Theorem 2, we find for An = t x X e , and Sn = — 2 rz that the 

termination condition 5 n = A n s n _i — A n _is n = yields 

1 



(ax 2 + 6x + cp 1 n + (/C ~ 1)a) ~ ^ n = 1 ' 2 ''-- ( 13 ) 

which yields for 5 n = that /i n = n(d + (n — l)a) as required. For n = 0, 1, 2, . . . i.e. 
/in = 0, fii = d, Hi — 2(d + a), . . ., we obtian by 

yn = exp(- Jj^dt), n = 0,l,2,..., (14) 
the polynomial solutions just mentioned. □ 

In Table I, we summarize the well-known differential equations which have 
polynomial solutions (as eigenfunctions). In each case, we give the explicit criterion, 
8 n — 0, of Theorem 2. 
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Table I: Application of AIM to classical differential equations. For each differential 
equation which give the condition under which it have polynomial solutions. 



3.1. Some remarks on Table I 

1 This differential equation is a generalization of the Cauchy-Euler linear equation 

x 2 y" + axy' + (3y = 0. (15) 

It is possible, however, to apply AIM to the differential equation (fToT) . The 
termination condition yields in this case 

( — 1)^+1 n 

5 n= 2w+2 \{{l3 + i{l-i + a)) = or = n{n - 1 - a) (16) 

x »=l 

while the corresponding polynomials, as given by (fTlj) . are y = 1> 2/i = ^> 
y 2 = x 2 , . . . ,y n = x n . It is clear that these polynomials cannot form an orthogonal- 
polynomial sequence [TU] . 

26 This differential equation can be regarded as a generalization of the well-known 
Hermite differential equation 2 ". It is an elementary example of differential equation 



Criterion for polynomial solutions 



6 



la 



with non-rational coefficients (i.e. with s and Ao non-rational) which has 
nonconstant polynomial solutions for 0. 

3 This is known as the confluent hypergeometric differential equation. It is also known 
as Kummer's differential equation or Pochhammer-Barnes equation [TT?] . 

6 This differential equation is known as Chebyshev's differential equation of the first 
kind and Chebyshev's differential equation of the second kind, respectively. It is 
interesting to note that these differential equations are special cases of 

{1 - x 2 )y" - axy' + fiy = 0. (17) 

If we apply AIM directly to (|17jl. we have by means of the termination condition 
(fTTl) that 

1 n 

Y[(i(i + a- 1) 



(- -')"•', 

thus, for S n = 0, we must have /i n = n(n + a — 1). The corresponding polynomial 
solutions, for n = 0,1,2,..., are y = l,yi = x,y 2 = (a + l)x 2 — l,y 3 = 
(a + l)x 3 — 3x, . . . , and in general 

CL 1 — X 

Vn = 2^i(-n, n + a-1,-, — ^— ) 
up to a constant. Here, 2 -Pi, Gauss' hypergeometric function, is defined by 



2 Fi(-n, b; c; x) = ^ ' 7 {" : . > 



c) fr A;! 



where the Pochhammer symbol (a)^ defined by 

Y(a + fc) 

(a) = 1, (a)fc = a(a + l)(a + 2) . . . (a + fe - 1) = —. . — . 

T(a) 

5a,fe rj^ e p iy nom i a ] solution of these differential equations were studied by Krall and 
Frink 20] . The corresponding (Bessel) polynomial solutions are orthogonal in the 
quasi-definite sense [TUj . 

In Table II we find the corresponding polynomial solutions for each differential equation 
mentioned in Table I. As an elementary application to quantum mechanics, we consider 
the Schrodinger equation 

+ ( A , 7(7 + 1)' 
dr 2 \ r r 2 

Writing ip(r) = r 1+1 e~ ar y(r), we easily find that y(r) must satisfy, for E = —a 2 , the 
confluent hypergeometric differential equation 

y"(r) =2[a- ^-jy'ir) + [ J -)y{r). 

The termination condition, mentioned in Table I, then yields E = —a 2 = — ^ n _^ +1 y , 
the eigenvalues of Schrodinger's equation for the Kratzer potential. Furthermore, the 
corresponding (un-normalized) eigenfunctions are given, by means of Table II, as 

Mr) = (-l) n ^ +1 e-^ r (2 7 + 2) nl F x (-n; 2 7 + 2; 2^Er). 
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solutions for each differential equation 



DE 


y n , n = 0, 1,2, . . . 


Cauchy-Euler 


yo = i 

y 1= x-b 

y 2 = (a- l)(a - 2)x 2 + 2{a - l)(2a - ab)x 
+ a 2 b 2 - a(2b + a) a + 2a 2 


Hermite 


yo(x) = l 
yi(x) = x 

y 2 (x) = 2x 2 -l 

y2n(x) = (-1)"2" (l/2) n l F 1 (-n; 1/2; x 2 ) , 
y2n+i(x) = (-l) ra 2"(3/2) n x iF!(-n; 3/2; x 2 ) 


Hermite 


yo( x ) = i 

2/i(x) = ax + 6 

2/2 (#) — + o) 2 — a 

y 2n (x) = (-l)"(2ar(l/2)„ ^ (-n; a/2; (ax + 6) 2 /2) , 

y 2n+1 (x) = (-l)"(2a) n (3/2) n (ax + 6) ,F, {-n; 3a/2; (ax + 6) 2 /2) . 


Laguerre 


2/0 = 1 

y 1 = x - 1 

y 2 = ^ 2 — 4x + 2 
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Hypergeometric 


2/o = 1 
yi = x + c 

y 2 = 2x 2 + 4(c + l)x + c(c + 1) 
2/n = (c)„ 2F (-ra, -n; c, x) 
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Generalized Bessel 
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j/i (x) = ax + b 

y 2 (x) = (a + l)(a + 2)x 2 + 26(a + l)x + 6 2 
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y n (x) = b n 2 F (-n, n + a-1;-; -x/b) 



Criterion for polynomial solutions 



9 



3.2. The case of A = 

In the early development of the asymptotic iteration method ^l] , one get the impression 
that the method is not applicable in the case of Ao = 0. This impression naturally arises 
because of the condition = n = 0, 1, 2, . . .. If Ao = 0, however, we may have 

using (JHJ) and (jBJ) that = — - — for which the corresponding asymptotic 

b n — l 

condition now reads 

hi = hzl = ai n = l,2,... (19) 

This leads to the essentially equivalent termination condition (|TT|) 

0~n = KSn-1 - K-lSn = 0, 71=1,2,... 

A simple example which show the use of AIM in case of Ao = is the differential equation 
x 2 y" — 2y = 0. Direct use of AIM implies that 62 = and a polynomial solution by 
mean of ()14j) is y 



x 



2 



4. Application to generalized Hermite, Laguerre, Legendre and Chebyshev 
differential equations 

Theorem 4: For N a positive integer and a,b ^ 0, the second-order linear differential 
equation (known as the generalized Laguerre differential equation) 

u " = (ax N - ^)u' - acx N ~ 1 u, (20) 

has a polynomial solution if 

c = n(N + l), 71 = 0,1,2,.... (21) 

The corresponding polynomial solutions are 

/ & + jV\ h + N ax N+1 



Proof: For N = 1, the termination condition (jTTj) yields c = 2n, n — 0, 1, 2, . . . while 
(JT3J) implies 

1, if n = (or c = 0) 

1 + b — ax 2 , if n = 1 (or c = 2) 

3 + 46 + 6 2 - 2a(3 + o)x 2 + aV, if n = 2 (or c = 4) 



2 n (^r^ ^(-n;^^;^), for n = 0,l,2,... (or c = 2n) 



2 Jn ±y ' 2 ' 2 



Criterion for polynomial solutions ... 10 

For N = 2, the termination condition (|11JI yields c = 3n, n = 0,1,2,... while (|14j) 
implies 



u n (x 



1, 

2 + 6 - ax 3 , 

10 + 76 + b 2 - 2a(5 + 6)x 3 + aV, 



'& + 2\ 6 + 2 ax 3 . 

3 [ —b Fi ^—'^ 



if n = (or c = 0) 
if n = 1 (or c = 3) 
if n = 2 (or c = 6) 

for n = 0, 1, 2, . . . (or c = 3n) 



Similarly, for iV = 3, the termination condition (fTTj) yields c = 4n, n = 0, 1, 2, . . . while 
(fT4]l implies 



1, 

3 + b — ax 4 , 

21 + 106 + 6 2 - 2a(7 + 6)x 4 + aV 



,., ,'b + 3\ ^ , 6 + 3 ax 4 , 

4| — J „ lFl( -" ; — ; - ; 
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if 
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= 1 
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= 4) 


if 
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= 2 


(or 
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= 8) 



for n = 0, 1, 2, . . . (or c = 4n) 



Similar expressions can be obtain for N = 4,5, . . .. These results can be generalized by 

Theorem 5: If b = 0, the second- order linear differential equation A2U\) . known as 
generalized Hermite differential equation, has a polynomial solution if 



c = n(N+l), n = 0,1,2, 



or 



c = n(N + l) + l n = 0,l,2,... 
In case of c = n(N + 1), the polynomial solutions are (for n = 0, 1, 2, . . .) 



Un (x) = (-ir(i\r+ir( T 



N 



M-n; 



N ax 



N+l 



+ Nj n 1 + N' 1 + N' 

In case of c = n(N + 1) + 1, the polynomial solutions are (for n — 0, 1, 2, . . .) 

u n {x) = (—1) (N + 1) x iFi(-n; ; ). 

nK J y 1 y 1 \\ + N) n 1 U 'l + iV'l + A^ ; 

Proof: Similarly to the proof of Theorem 4, the conditions (|23j) and ()24|) follow directly 
by means of the termination condition (jllj) . with Ao = ax N and s = — acx^^ 1 . Eqs 
(|23|) and © follow from JTJ). □ 



(23) 
(24) 

(25) 
(26) 



Theorem 6: For N a positive integer, the differential equation 



ax 



N 



wx 



N-l 
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has polynomial solutions 



u n (x) 



(N+l)- n \N+l) 





+ n; 



b + N 
1 + N 




(28) 



w = n (N + l)(s(6 - 1 + n(N + 1)) + a), n = 0, 1, 2, . . . 



(29) 



where 2 -Pi 2s Gauss's hypergeometric function i!8\) . If b = 0, s = 1, i/ien £/ie cases 
of a = 2 and a = 1 corresponding to differential equations known as the generalized 
Legendre and Chebyshev differential equations, respectively. 

Proof: Using AIM, condition (J2U|) for polynomial solutions follows by means of the 
termination condition 8 n = in a similar fashion to the proof of Theorem 4. Equation 
(j2Hj) then follows by means of (fT^j) as generalization of the polynomial solutions for each 



5. Conclusion 

We have presented a simple criterion for the existence of polynomial solutions of second- 
order linear differential equations. Many of the classical differential equations that 
appear in mathematical physics can be analysed with this theory. Apart from its 
theoretical interest, the criterion can be used in a practical way to look for and to 
obtain polynomial solutions to eigenvalue problems of Schrodinger-type JIJ-JH], and 
similarly for polynomial solutions of quasi-exact solvable models in quantum mechanics 
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